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1.0 INTRODUCTION 


This report summarizes work conducted over the past year designed 
to examine the mechanics of fatigue crack growth under constant-amplitude 
and variable-amplitude loading. In the latter case, critical loading 
histories involving relatively simple overload and overload/underload 
cycles were studied in order to provide a basic understanding of the under- 
lying physical processes controlling crack growth. The material used 
for this study was 7091-T7E69, a powder metallurgy aluminum alloy. Exten- 
sive measurements of local crack-tip parameters have been made at various 
times before, during, and after the overloads--these include crack-tip 
opening loads and displacements as well as crack- tip strain fields. The 
latter were used, in combination with the materials cyclic and monotonic 
stress-strain properties, to compute crack-tip residual stresses. The 
experimental results were also compared with analytical predictions obtained 
using the FAST-2 computer code. A study of the sensitivity of the analyti- 
cal model to constant-amplitude fatigue crack growth rate properties and 
to through- thickness constraint was also conducted. 

The report is organized into four main sections which are distinct 
but interrelated. 

Section 2.0 presents crack closure measurements as a function of 
the applied AK and K max , as well as a comparison of crack closure measure- 
ments obtained local to the crack tip using an optical stereoimaging tech- 
nique and remote from the crack tip using an elastic compliance technique. 

Section 3.0 contains a detailed determination of the residual stress- 
strain fields at crack tips resulting from overloads and overload/under- 
loads. These results are used to establish a correlation between crack-tip 
strain and crack growth rate and to examine the issue of similitude by 
comparing stress-strain fields obtained before, during, and after the 
overloads . 

Section 4.0 presents the sensitivity study of the FAST-2 computer 
model and compares analytical predictions with experimental results. 
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Section 5,0 provides a general discussion of results from the entire 
study, and Section 6.0 gives the resulting conclusions and recommendations. 

Two appendices are provided for supplemental information. Appendix 
A gives the definitions of the terms used to characterize the overload 
and overload/underload events, as well as resulting crack growth retardation 
and delay phenomena. Appendix B provides a supplement to Section 5.0 
by presenting additional details on the comparison of analytical and experi- 
mental results, as well as predictions of the crack-tip residual stresses 
from the FAST-2 analyses. 

The work summarized in this report is an extension of that conducted 
under a previous NASA program [1]. 
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2.0 CRACK CLOSURE MEASUREMENTS FOR CONSTANT-AMPLITUDE FATIGUE CRACK GROWTH 


Extensive measurements have been made of the load at which a fatigue 
crack opens, for the powder metallurgy aluminum alloy 7091. These measure- 
ments were made on single-edge notched specimens, Figure 1, with the stereo- 
imaging technique, using both an optical imaging system and a scanning 
electron microscope (SEM). Also, for the 2024 center-cracked panels sup- 
plied by NASA shown in Figure 1, measurements of crack opening load were 
made both locally and remotely. Locally, they were made using optical 
microscopy and stereoimaging; remotely, the compliance technique was used, 
and potential drop was attempted. Results from the various methods are 
compared. 

2. 1 Local Crack-Tip Opening Load as a Function of AK and K^ ^y 

Crack opening load is determined using the stereoimaging technique 
by comparison in a stereoviewer of photographs made at successively-increas- 
ing loads with a photograph made at the minimum cyclic load. The viewing 
axis is that of Mode I; i.e., the eye axis is the same as the loading 
axis. Crack opening load in Mode II can also be determined by imaging 

the crack-tip region with the eye axis orthogonal to the loading axis. 
The latter procedure was used in a limited number of cases. 

When observed in this way, the crack is seen to "peel open" with 

increasing load until it is open down to the crack tip. The load required 
to obtain complete opening of this crack to its tip in Mode I is termed 
P Q I. The accuracy of this method of determining opening load is greater 

than its reproducibility from cycle to cycle. That is, if opening loads 
are measured carefully each 10-100 cycles at constant AK, variations in 
magnitude result; this observation is consistent with the previous findings 
of periodic variation in crack-tip opening displacements and strains as 

the crack propagates. 

Typical findings of the variation in the length of the crack which 
remains closed with increasing applied load are shown in Figure 2. The 
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(a) 7091 aluminum alloy (Alcoa) 



FIGURE 1. SPECIMEN DESIGNS USED. (a) Single edge notched design used for 
for both constant amplitude and spectrum loading using both 
scanning electron and optical microscopy, (b) Center notched 
panel for crack opening load measurement under constant ampli- 
tude loading as determined by stereoimaging and compliance. 
Dimensions in mm. 




Applied Load Fraction, p /Pmax 


7091 

R = 0.1 


Arrows indicate Pp^Pmax 


AK = 6 MN/m3/y° 


AK = 10 


Distance behind crack tip, -y (um) 


FIGURE 2. CRACK OPENING LOCATION RELATIVE TO THE CRACK TIP AS LOAD P IS 
INCREASED TO Pq 1 , THE CRACK OPENING LOAD IN MODE I. Note the 
large change in load required to open the crack the last few 
micrometers at low AK. 
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spatial resolution of the measurement technique is clearly an important 
factor in determining the magnitude of the opening load because of the 
large change in values close to the crack tip. This is especially true 
for values of AK near the threshold for fatigue crack growth. 

The results of numerous measurements obtained over a wide range 
of applied AK and R are presented in Figure 3. Values at low-R and low-AK 
were made in the SEM under high-resolution conditions, while values at 
high-R were made optically at lower resolution. A comparison between 
optical measurements and SEM measurements at AK = 6 MN/m 3 / 2 and R = 0.1 
showed that optical microscopy resulted in AK e ff values 15% larger than 
those determined in the SEM, whereas at AK = 10 MN/m 3 / 2 , the results were 
coincident. At higher R, optical and SEM-derived results were also in 
agreement, at least down to AK = 6 MN/m 3 / 2 . 

The values of AK e ff / AK in Figure 3 are calculated from measured 

values of P 0 /P max us i n E the equation 

AK eff /AK = (1 - P 0 /P max )/ ( 1 - R) (1) 

The line shown on the figure is not fit statistically because of the large 
scatter in the data and because some data are better than others due to 
the learning process which took place as more data were obtained. Neverthe- 
less, the line shown is a fair representation of the results. The intercept 
of the line at AK e ff/AK = 0 is defined as the threshold for fatigue crack 
growth, and is shown as 5 MN/m 3 / 2 for this case. Values measured from 
crack growth experiments are known to be lower than this, which complicates 
interpretation of this concept. The reason that the threshold is lower 
than 5 MN/m 3 / 2 is believed to be related to the fact that the Mode II 

opening load is lower than the Mode I opening load at low AK. Thus, Mode 
I no longer controls the value of AK; rather it is Mode II which sets 
the lower limit of fatigue crack growth. 

The results of Figure 3 may be used to reinterpret the comparison 

between measured and computed values of AK e ff/AK and R. The comparison, 
shown in Figure 4, indicates that a family of lines should be used to 
show different values of AK in this representation. Results of FAST-2 
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FIGURE 3. RELATION BETWEEN NORMALIZED EFFECTIVE STRESS INTENSITY FACTOR 
AND Kniax, AS DETERMINED FROM LOCAL MEASUREMENTS OF OPENING 
LOAD IN MODE I. 



AK eff /AK 



0 0.2 0.4 0.6 0.8 1.0 


^min/^max 


FIGURE 4. COMPARISON OF MEASURED AND PREDICTED NORMALIZED EFFECTIVE CRACK-TIP 
STRESS INTENSITY RANGES (AK eff /AK) AS A FUNCTION OF LOAD RATIO (R) 
FOR CONSTANT-AMPLITUDE FATIGUE CRACK GROWTH. 
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were obtained at AK = 7 MN/m^/Z and are in good agreement at intermediate 
R values, but diverge at both high and low values of R. At high-R, the 
measurements are extrapolated, so the uncertainty between measurement 
and computation is not resolved; however, at low-R, the measurements are 
very reliable and the divergence between computation and experiment is 
meaningful . 

2.2 Comparison of Compliance and Stereoimaging Measurements of Crack 

Opening Load 

Panels of 2024-T3 aluminum alloy (2.3 mm thick) supplied by NASA 
were used for these measurements. The specimen design is shown in Figure 
1(b). Both potential drop and compliance techniques were used for remote 
measurements. An optical microscope, together with the stereoimaging 
technique, was used for the local measurements. 

For the potential drop measurements, current leads were attached 
by spot welding near the loading pin holes and the potential leads were 
attached at the small (4.8 mm diameter) holes on either side of the center 
notch. Autographic recordings of potential versus applied load provided 
no indication of crack closure. It is presumed that oxidation of the 
crack surfaces due to moisture in the air was the apparent cause of the 
failure to obtain crack opening measurements by potential drop, and efforts 
to make meaningful measurements with this technique were abandoned. In 
order to have been successful, it appears that it would have been necessary 
to conduct these experiments in a vacuum or inert gas environment. 

All compliance measurements were made with the compliance gage 
located 2 mm behind the crack tip. Selected measurements were also made 
with the compliance gage located 8 to 10 mm behind the crack tip. However, 
no measurable difference in crack opening load was obtained for these 
two locations. For the compliance gage used, the gage length of the dis- 
placement measurement was 5 mm. 

Fatigue cracks were extended under constant-amplitude loading at 
R = 0.05 and P max = 1053 lb, and measurements were made at relative crack 
lengths (a/W) of 0.26. 

Results of crack opening load measurements by compliance are compared 
with those by stereoimaging of optical photographs in Figure 5. Both 


Comp. S. I. 

(a) 


Comp. S. I. 

(b> 


FIGURE 5. COMPARISON OF THE OPENING LOAD TO MAXIMUM LOAD RATIO AS 
DETERMINED REMOTELY (COMPLIANCE) AND LOCALLY (OPTICAL 
MICROSCOPY WITH STEREOIMAGING) FOR (a) 2024 PANELS SUP- 
PLIED BY NASA, AND (b) 7091. (a) For 2024, both measure- 

ments were made at SwRI on the same specimen, while (b) 
for 7091, local measurements made at SwRI are compared 
with compliance measurements reported by McEvily. 
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compliance and stereoimaging measurements were made on the 2024 specimens 
at SwRI, while for the 7091 measurements, the compliance values were derived 
from work reported by McEvily [2]. As the figure illustrates, stereoimaging 
measurements in each case give higher values of opening load than do the 
compliance measurements, a finding consistent with the results of Figure 
2 . 
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3.0 CRACK-TIP STRAIN AND STRESS FIELDS DURING CONSTANT-AMPLITUDE 
AND VARIABLE-AMPLITUDE FATIGUE CRACK GROWTH 

Three overload sequences have been analyzed in great detail, with 
one of the overloads being followed by an underload. In each case, measure- 
ments were made before, during, and after the overload event; thus informa- 
tion was obtained for both constant-amplitude and variable-amplitude fatigue 
crack growth. The history of each load excursion analyzed is shown in 
Figure 6, and the relevant loading parameters for each case are listed. 
A comparison of the crack growth rates subsequent to each of these load 
excursions is shown in Figure 7, along with determinations of effective 
stress intensity factors. 

The methodology used for analysis of each load excursion will be 
detailed for one of the overloads in the next section. Summary information 
for all cases examined will then be presented. 

3.1 Analysis of Experiment 1: OLR* = 2.15, ULR* = 1, R = 0.16, AK = 

6.2 mh? n 

In preparation for the load excursion event, a photograph was taken 
under minimum load at a location several hundred micrometers ahead of 
the crack tip. This is labeled as Photograph No. 1 on Figure 8. The 
specimen was then cycled until the crack had grown into the field of this 
photograph, and a matching picture (No. 2) was taken at minimum load. 
Next, Photograph No. 3 was taken at the maximum load of the cyclic range. 
Comparison of Photographs 1 and 2 shows the effect of the crack being 
present, and comparison of Photographs 2 and 3 gives the cyclic crack-tip 
parameters just prior to the load excursion. The crack opening load was 
also measured on this last constant amplitude cycle. 

The overload cycle was then applied, with Photograph No. 4 taken 
at the maximum load, and No. 5 taken at the minimum load, as shown in 
Figure 8. Comparison of Photographs 2 and 4 shows the effect of the loading 
portion of the overload, while the reversibility due to the subsequent 
unloading part of the cycle is shown by comparison of Photographs 4 and 
5. 
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Note: For each sequence stereoimaging data were taken 

for both loading and unloading of overload cycle 
plus two or three additional times during the 
subsequent delay period. 


FIGURE 6. DEFINITION OF LOADING PARAMETERS AND SUMMARY OF SEQUENCES 
FOR WHICH RESIDUAL STRESSES WERE DETERMINED. 
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iA roll owing Over! old, mn 

(a) Measured crack growth rate response 
and AKgff following a single over- 
load at R =0.16, OLR* = 2.15. 



la Following Overload. m 

(b) Measured crack growth rate response 
and AKgff following a single over- 
load at R = 0.22, OLR* = 2.85. 



11 Following Overload, m 

(c) Measured crack growth rate response and AKeff 
for a single overload followed by an underload 
to near zero (R = 0.5, OLR* = 3, OLR* = 2). 


FIGURE 7. CRACK GROWTH RATE AND EFFECTIVE STRESS INTENSITY FACTOR HISTORIES OF THE 
THREE LOAD EXCURSIONS ANALYZED USING STEREOIMAGING. Cases (a) and (c) 
are similar to those for which results were obtained using FAST-2* 


AK efr /aK 
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Load, N 



113B 114A,B 118B 119A,B 122 123A,B 123C,D 


AK = 6.9 AK = 6 AK = 7.2 

R = .16 R = .22 R = .5 

OLR* = 2.15 OLR* = 2.85 OLR* = 3.0 

ULR* = 1 ULR* = 1.0 ULR* = 2.0 

FIGURE 8, LOADING HISTORY FOR THE THREE LOAD EXCURSIONS ANALYZED, The 
circled numbers indicate the points where photographs were 
taken for subsequent analysis. Below each part of the figure 
is information describing the loading history. 



17 


Data sets are labeled as indicated in Table 1 and partially shown 
in Figure 8. Examination of these data were made using the stereoimaging 
technique [3], and measurement of the displacements was accomplished using 
the methods of photogrammetry . Strains were computed from the gradients 
of the measured displacements [4]. 

3.1.1 Results Before the Overload 

The net displacements caused by the growth of the crack (comparison 
of Photographs 1 and 2) are shown in Figure 9. The displacements are 
seen to be small and are local to the crack region, as anticipated. Most 
of the displacements occur within a few micrometers of the crack. Net 
displacement in the loading direction across the region marked Dr on Figure 
9 is 0.12 urn. Strains have been computed for this case and are shown 
in Figure 10 as Mohrs circles. With a few exceptions, the strains are 
less than a few percent. 

For constant-amplitude cyclic loading prior to the overload, 
displacements (Figure 11) are similar to previous work with this material 
[5]. The distribution of maximum shear strains in the near crack tip 
vicinity is shown in Figure 12. Note the asymmetric and biaxial nature 
of the displacements about the crack, but' the more symmetric distribution 
of shear strains. These findings are similar to previous results for 
7075-T651 and T1-6A1-4V. 

3.1.2 Results for the Overload Cycle 

Displacements caused by the loading portion of the overload cycle 
are shown in Figure 13, and those due to unloading are shown in Figure 
14. Greater symmetry is evident in the overload cycle than is found in 
the cyclic response prior to the overload, and considerable reversibility 
in displacements is evident upon unloading. 

Strains computed from the displacements shown in Figures 13 and 
14 have the same general distribution as the cyclic strains before the 
overload, but peaked much more near the crack tip. Further analysis of 
these strains will be found in the next section. 
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TABLE 1 

LABELING OF DATA SETS ANALYZED 


Suffix 

Preoverload 

A 

B 


Condition 


no crack vs unloaded crack 
photograph 1 vs 2 

unloaded crack vs loaded crack 
photograph 2 vs 3 


Overload Cycle 
A 

B 


loaded crack vs maximum load 
photograph 3 vs 4 

crack at maximum load vs unloaded crack 
photograph 4 vs 5 


Underload Cycle 
C 


crack unloaded after overload vs maximum 

underload 

photograph 5 vs 6 

crack at maximum underload vs unloaded 
crack 

photograph 6 vs 7 



FIGURE 9 
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DISPLACEMENTS CAUSED BY THE PRESENCE OF THE CRACK. These were 
determined by comparison of photographs #1 and #2 of Fig. 8. 
Note that displacements were measured each two micrometers 
(scale on the left), and that the scale of displacements is 
greatly magnified (scale on the right). The dot marks the 
original position of the material point, and the tip of the 
line is the point to which the material has been displaced by 
passage of the crack. 
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2 MICRONS 1 0 7. STRAIN 

FIGURE 10. MOHRS CIRCLES OF STRAIN AS DERIVED FROM THE DISPLACEMENTS SHOWN 

IN FIGURE 9. Maximum shear strain is the diameter of the circle, 
while the mean strain is indicated by a numeral. 
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FIGURE 11. DISPLACEMENTS FOR THE CYCLIC CRACK JUST PRIOR TO THE OVERLOAD 
CYCLE. Note that the displacement magnitude is larger than 
the actual spatial scale. 
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FIGURE 12. DISTRIBUTION OF THE MAXIMUM SHEAR STRAIN FOR THE CYCLIC CRACK JUST PRIOR TO THE OVERLOAD 
CYCLE. The crack is shown schematically on the plane of zero strain. 
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FIGURE 13. DISPLACEMENTS CAUSED BY THE CRACK BEING LOADED TO MAXIMUM 
VALUE DURING THE OVERLOAD CYCLE. These displacements were 
determined from comparison of photographs #2 and #4 of 
Figure 8. 
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20-00 MICRONS 


6-00 MICRON 
DISPLACEMENT 


FIGURE 14. DISPLACEMENTS CAUSED BY THE UNLOADING PORTION OF THE OVERLOAD 
CYCLE. Comparison with Figure 13 indicates considerable re- 
verse material flow. 
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3.1.3 Computation of Stresses 

One of the major efforts of this analysis has been to compute the 
distribution of residual stresses which result from the load excursion. 
Stresses were computed from strains at each step in the applied loading 
sequence, using a procedure previously derived [6]. This computation 
requires knowledge of the appropriate stress-strain curves. Consequently, 
the monotonic and cyclic stress-strain curves were measured for this mate- 
rial and were reported previously [1]. The cyclic and monotonic curves 
are nearly equal and there is essentially no difference between the tension 
and compressive monotonic curves. The cyclic hysteresis loops were nearly 
perfect in symmetry, and stabilized very rapidly after a change in strain 
range . 

To obtain the residual stress values, coordinate values of each 
stress, together with the shear stress, were computed at each node around 
the crack tip [6], for each step in the loading sequence, and then added 
together to produce a net residual stress, which occurred due to the nonuni- 
formity and nonreversibility of deformation within the plastic zone. 

For each load excursion, the result of the stress computation was 
six components of stress at each point of the analysis. Two of these 
stresses were chosen for detailed analysis: the stress in the loading 
axis (o x ) and the total effective stress (a e ff) were chosen because the 
loading axis stress could be expected to exert a direct clamping influence 
on the crack, and the effective stress was examined because it considers 
both the normal and shear stresses in its computation and represents a 
stress equivalent to what might be applied in a tensile test. 

3.1.4 Determination of Residual Stress 

Mostly positive residual stresses in the loading axis resulted, 
in each case analyzed, from summing the stresses on the loading and unload- 
ing portions of the overload cycle, as would be expected, because the 
net deformation caused by the overload was generally tensile, and no compen- 
sation was made for the fact that the overload cycle began at the minimum 
in the hysteresis loop, which is at a negative stress. The situation 
described above is depicted in Figure 15(a), while Figure 15(b) shows 



FIGURE 15. DISTRIBUTION OF THE STRESSES RESULTING FROM THE OVERLOAD. Th. 

stresses resulting from the loading and unloading portions of 
the cycle were added to produce the stresses shown. (a) is 
the distribution of stresses in x, the axis of loading, and 
(b) the distribution of stresses in y, the direction of crack 
growth. 
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the stress perpendicular to the loading axis, and Figures 15(c) and (d) 
show the maximum principal stress and the effective stress, respectively. 
In order to adjust the stresses to their correct value along the stress 
axis, it was thus necessary to subtract the value of the negative residual 
stress caused by the previous cyclic loading history of the crack. 

The stresses in the specimen at the minimum cyclic load and caused 
by the constant-amplitude cyclic loading history can be derived from using 
the data of Figure 11, which resulted in analysis of Photographs 2 and 

3 (data sets "A”). First the stresses were computed from the strains, 
but because of the localized nature of the cyclic loading, these data 
were taken from photographs at higher magnification than the photographs 
of the subsequent overload cycle. Therefore, in order to obtain common 
nodal points from which stresses could be added, a model of the stresses 
for the cyclic loading prior to the overload was constructed from the 

computer values. The modeling procedure consisted first of graphing the 

computed stresses. These were then made symmetric about the line of the 

crack because, on average, this is the behavior which has been experi- 
mentally found. Close to the crack tip, stresses computed from strains 
were used, but farther away from the crack, where plasticity had ceased 
to dominate, these stsresses were faired . into the stresses computed for 
an elastic crack. A comparison of the actual and modeled stresses is 
shown in Figure 16. 

This modeling procedure allowed the stresses to be computed at 
the same nodal points for which analyses were made after the load excursion. 
The stress field derived using this method is shown in Figure 17. The 
stresses shown are the change in stress between minimum and maximum cyclic 
load. 

The next step was to determine how much of the stress range computed 
was actually compressive when the load was at minimum. Opening load was 
used to make this estimate. The stress range was partitioned using the 
fraction of total K change which was effective, as computed from the opening 
load. The concept is shown in Figure 18. The fraction of stress computed 
to be negative from the cyclic loading analysis was then subtracted from 
the net stress resulting from the load excursion. The resulting stresses 



FIGURE 15 (CONTINUED). DISTRIBUTION OF THE STRESSES RESULTING FROM THE 

OVERLOAD. The stresses resulting from the load 
ing and unloading portions of the cycle were 
added to produce the stresses shown, (c) is the 
distribution of maximum principal stress, and 
(d) the distribution of the effective stress. 
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10 MICRONS 1100MPA STRESS 

(b) 


FIGURE 16. MOHRS CIRCLES OF STRESS FOR THE CYCLIC CRACK PRIOR TO THE 
OVERLOAD. (a) as determined directly from the measured 
strains, and (b) as determined from a model of the crack 
tip stress field. 
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10 MICRONS 


900 MPA STRESS 


FIGURE 17. CYCLIC STRESSES AS DETERMINED FROM THE MODEL AND EXTRAPOLATED 
TO COVER THE SAME FIELD AS THE OVERLOAD CYCLE ANALYSES. 
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Stress 




FIGURE 18. PARTITIONING OF THE CYCLIC STRESS RANGE PRIOR TO THE OVERLOAD. 

Partitioning is done by using the measured crack opening load, 
and is assumed to be equally distributed throughout the plastic 
zone. 
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in the load axis (a x ) and the effective stress are shown in Figures 19 
and 20. In Figure 21, a cross section through the residual stress field 
of Figure 19 is shown along the crack plane (x = 0), in order to make 
visualization of the effect clearer. 

3.1.5 Post-Overload Cycling 

After the load excursion, the cyclic loading range was returned 
to that previously applied, and crack growth was monitored as the crack 
progressed through the plastic zone of the load excursion. At two points 
along this growth path, additional analyses were made of the crack-tip 
strain field. One point was chosen to be as close as possible to the 
point of minimum crack growth rate. The opening load at each point was 
measured and the strains around the crack tip were determined. An analysis 
was made to determine if the crack tip strain fields were in similitude 
with the crack-tip field before the load excursion and whether the crack-tip 
strains correlated with crack growth rates. 

3.2 Comparative Analysis of the Load Excursions 

The analysis method outlined in the previous section was applied 
to all the loading variations shown in Figure 8. This section summarizes 
the findings of the other two load variations studied and compares the 
results . 

3.2.1 Overload Cycle Strains 

Maximum shear strains at the peak of the overload cycle are compared 
in Figure 22. It is apparent that strain magnitude increases with the 
magnitude of the stress intensity factor. Crack-tip strain fields have 
also been determined for a crack loaded monotonically in tension, and 
the crack-tip strains for the overload cycle and monotonic loading are 
compared in Figure 23, where it may be seen that the correlation between 
crack-tip strain and stress intensity factor for monotonic loading is 
different than for the overload cycle. 

The distribution of strains directly ahead of the crack has been 
examined for the three overload cycles and is shown in Figure 24. Strains 




FIGURE 19. THE RESIDUAL STRESS IN THE DIRECTION OF LOADING WHICH RESULTED FROM APPLICATION OF THE 
OVERLOAD. These stresses are the sum of those shown in Figure 15(a) and Figure 17, as 
partitioned using the concept of Figure 18. uj 







FIGURE 21. A CROSS SECTION THROUGH THE RESIDUAL STRESS FIELD OF FIGURE 19 SHOWING THE STRESSES 
DIRECTLY AHEAD OF THE CRACK TIP CAUSED BY THE APPLICATION OF THE OVERLOAD CYCLE. 
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FIGURE 22. COMPARISON OF THE MAXIMUM SHEAR STRAINS 
RESULTING FROM EACH OF THE OVERLOAD 
CYCLES ANALYZED. (a) OLR* = 2.15, (b) 
OLR* = 2.85, and (c) OLR* = 3.0. 
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FIGURE 23. COMPARISON OF THE EFFECTIVE STRAIN AT THE CRACK TIP CAUSED BY 
THE OVERLOAD CYCLE WITH THE STRAIN CAUSED BY MONOTONICALLY 
LOADING A CRACK FROM AK = 5 MN/m 3 / 2 . 




Normalized Strain 



FIGURE 24. OVERLOAD CYCLE EFFECTIVE STRAINS AHEAD OF THE CRACK TIP, AS NORMALIZED BY THE CRACK TIP STRAIN. 

The horizontal axis is the adjusted distance parameter (r+B). The line shown represents the 
equation shown in the figure. The fact that one equation fits all data is one indication that 
similitude exists between all three crack tips on the overload cycle. 
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from all three may be described by the following function: 

e' (r) = 1.54 ln(-^|) (2) 

where z'(r) = strain at the distance r, normalized by the crack-tip strain. 
The constant B = 0.33 and limits the strain to a finite value at the crack 
tip (r = 0). This result may be compared to a similar analysis of strains 
for monotonic loading, Figure 25. For this case, the function shown on 
the figure appears to describe the strains in the near-tip region, but 
not in the far-field region. 

A comparison of results from Figures 24 and 25 suggests that differ- 
ences occur in the strain distribution within the plastic zone for fatigue 
cracks which are overloaded and cracks which are monotonically loaded. 

3.2.2 Cyclic Loading Strains 

The similarity of strain distributions ahead of the crack tip before 
and after the load excursion is examined in Figure 26. Although there 
is considerable scatter in the data, it is concluded that similitude is 
demonstrated by the crack- tip strain fields, and that it is not altered 
by the load excursion. This conclusion is further strengthened by the 
correlation of crack-tip strains and opening displacements, Figure 27. 
The cyclic data fall mostly along the same line, but the monotonic tension 
and overload cycle data do not. 

3.2.3 Plastic Zone Size Ahead of Crack Tip 

The functions which describe the strain distributions ahead of 
the crack tip are summarized in Table 2. There is obviously not much 
commonality between the three loading conditions. Apparently, loading 
history is important in determining the strain distribution ahead of a 
crack tip. Use of the normalized strain functions to determine plastic 
zone size for an overload from the monotonic results would, therefore, 
lead to a large error. Plastic zone sizes before, during, and after the 
load excursions have been derived from the relevant strain distributions 
and are listed in Table 3. As expected, the plastic zone size for the 



Normalized Strain 



FIGURE 25. 


MONOTONIC LOADING NORMALIZED EFFECTIVE STRAIN VS ADJUSTED DISTANCE PARAMETER. The line shown 
fairly well represents the data, except at high values of K and at distances exceeding about 
70 micrometers. Comparison with Figure 24 indicates a clear difference. 
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FIGURE 26. COMPARISON OF THE EFFECTIVE STRAIN DISTRIBUTIONS BEFORE 
AND AFTER THE OVERLOAD CYCLE FOR THE THREE OVERLOADS 
ANALYZED. Similitude appears to be approximately pre- 
served in spite of the presence of the overload residual 
stress field. 
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FIGURE 27 . CRACK TIP OPENING DISPLACEMENT COMPARED WITH CRACK TIP 
EFFECTIVE STRAIN FOR CYCLIC, OVERLOAD CYCLE AND MONO- 
TONIC LOADING. Also shown is data previously determined 
for a similar alloy. Crack tip opening displacement is 
defined at a distance of one micrometer behind the crack 
tip. 
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TABLE 2 

EQUATIONS FOR NORMALIZED STRAIN AHEAD OF CRACK TIP 

r in pm 


Loading 


Cyclic 

Overload cycle 
Monotonic 


Equation 


e' 

e' 

e' 


.22 In 
154 In 
180 In 


( 


( 


( 


57.7 

r+.612 

219 


r+. 33 

64.7 
r+. 32 


) 


TABLE 3 

DISTANCE AHEAD OF CRACK TO ELASTIC STRAIN BOUNDARY 
(Plastic Zone Size) 


Experiment 1: 

.OLR* = 2.15, 

ULR* = 1, 

R = .16 


Set 

Crack 


AK 

r(ym) 

113B 

Before 


6.2 

42.8 

114A 

OL 


14.8 

164 

115 



6.2 

29.3 

116 



6.2 

30.0 

Experiment 2 : 

OLR* =2.85 

, ULR* = 1, 

CM 

CM 

• 

II 

& 


118B 

Before 


6 

55.5 

119A 

OL 


17.1 

172 

120 



6 

23.7 

121 



6 

15.5 

Experiment 3: 

OLR* = 3, ULR* = 2, R 

= .5 


122B 

Before 


7.2 

46.6 

123A 

OL 


21.7 

190 

124 



7.2 

35.8 

125 



7.2 

41.9 
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overload cycle is larger than that for the pre-overload cycle, and the 
latter is larger than for cracks in the overload-affected period. 

3.2.4 Correlation Between Crack Growth Rates and Crack-Tip Strains 

The relationship between total crack-tip strain range and crack 
growth rate has been examined, Figure 28, for cracks growing before and 
after the load excursion. With one exception, which was explained, the 

correlation was excellent. Note that the results also correlate well 
with the severity of the load excursion. 

This result is another indication that similitude is satisfied 
after application of the load excursion. Similar conclusions were drawn 
from the results of Figures 26 and 27. Since the crack is responding 
in a predictable way, this indicates again that it is the driving force 
which is unknown subsequent to the load excursion. 

3.2.5 Residual Stresses from the Load Excursion 

Stress changes computed from the strain changes measured at each 
step in the load excursion are shown in Figures 29, 30, and 31, where, 

as indicated in Table 1, the "A" refers to the overload-loading part of 
the cycle and "B" refers to the overload-unloading subcycle. For the 

overload/underload cycle, Figure 31, "C" is the underload-unloading part 

of the cycle and "D" refers to the underload-reloading part. The results 
are about as might be expected from a knowledge of the other load excur- 
sions, except for the underload-unloading part of the cycle, Figure 31(d), 
where the stresses are more inhomogeneous than in the rest of the load 
excursion. 

The summed stresses, adjusted to a reference of zero stress, as 
outlined in the first section, are shown in Figures 32, 33, and 34. These 
are the residual stresses in the axis of the loading (o x ). Figures 32 

and 33 should be compared with Figure 19. Note that the overload from 

R = 0.5 is treated both before and after the subsequent underload. The 
three overloads yield similar results, although Figure 19 exhibits more 
inhomogeneity than the other two. A cross section through the residual 
stress field is shown along the crack direction in Figure 35, which should 
be compared with Figure 21. 


t/ Total crack tip effective strain 
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FIGURE 29. (a) STRESSES ON THE LOADING PORTION OF THE OVERLOAD CYCLE AND 

(b) ON THE UNLOADING PORTION OF THE OVERLOAD CYCLE IN THE 
AXIS OF THE LOAD APPLICATION, OLR* = 2.15. 
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FIGURE 30. (a) STRESSES ON THE LOADING PORTION OF THE OVERLOAD CYCLE AND 

(b) ON THE UNLOADING PORTION OF THE OVERLOAD CYCLE IN THE 
AXIS OF THE LOAD APPLICATION, OLR* = 2.85. 
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FIGURE 31. (a) STRESSES ON THE LOADING PORTION OF THE OVERLOAD CYCLE, AND 

(b) ON THE UNLOADING PORTION OF THE OVERLOAD CYCLE. All 
stresses are in the axis of load application. OLR* = 3.0. 
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FIGURE 31 (CONTINUED). (c) STRESSES ON THE RELOADING PORTION OF THE 

UNDERLOAD CYCLE, AND (d) ON THE RELOADING 
PORTION OF THE UNDERLOAD CYCLE. All stresses 
are in the axis of load application. 

ULR* = 2.0. 


SIGMA 


FIGURE 32. RESIDUAL STRESS FIELD IN THE DIRECTION OF LOADING CAUSED BY THE OLR* = 2.85. 
figure should be compared to Figure 19 (OLR* = 2.15). 


This 




2 °Oo. 

1 50 0 . 


1 

<0 



FIGURE 33. RESIDUAL STRESS FIELD IN THE DIRECTION OF LOADING CAUSED BY THE OLR* = 3.0. This figure 
should be compared with Figures 19 and 32. Note that the crack was not subsequently 
grown through this condition, but was followed by an underload. 



FIGURE 35. RESIDUAL STRESS DISTRIBUTION WITHIN HALF OF THE CRACK TIP PLASTIC ZONE FOR OLR* =3.0 AND 
R = 0.5. The section along X = 0 shows the profile of residual stresses normal to the 
crack whose tip is at Y = 0. Both X and Y are in pm and SIGMA X is in MPa. 
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The residual effective stresses appear similar to the residual 
stresses in the loading axis, Figures 36, 37, and 38. These should also 
be compared to Figure 20. 

The residual stresses of Figure 34 or 37, which are a result of 
the overload/underload, show both tensile and compressive stresses, result- 
ing in a very inhomogeneous residual stress field, but one which has a 
nearly zero net change in stress averaged over the plastic zone. This 
small magnitude of residual stress is reflected in the small change in 
crack growth rates as the crack moves through the residual plastic zone. 
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FIGURE 38. RESIDUAL EFFECTIVE STRESS FIELD FOR OLR* = 3.0 AND ULR* = 2.0. Compare 
with Figures 34 and 36. 
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4.0 ANALYSIS OF SELECTED EXPERIMENTS WITH THE PLASTIC WAKE MODEL (FAST-2) 

Three experiments involving either a single overload or an overload/ 
underload sequence were selected for detailed analysis with Newman's plastic 
wake closure model (FAST-2) [7]. Two of these experiments (shown in Figure 
7), designated Experiments 1 and 3, have been previously discussed in 
the detailed plasticity analysis of Section 3.0. One other experiment, 
designated Experiment 4, was also modeled. Table 4 shows the data charac- 
terizing these experiments. 

TABLE 4 

LOADING VARIABLES FOR EXPERIMENTS ANALYZED WITH FAST-2 

AK = 7 MPa4i 

Experiment 

No. OLR* ULR* R 

1 2.15 1.0 0.1 

3 3.0 2.0 0.5 

4 3.0 1.0 0.5 

The sensitivity of the analysis to variations in the degree of 
constraint and to variations in constant-amplitude fatigue crack growth 
properties were examined. The model was evaluated by comparing analytical 
predictions with experimental measurements. Specifically, the following 
quantities were compared for times before, during, and after the load 
excursions s 

• Crack opening stress, expressed as the effective 
load ratio (R eff = P 0 /P max ) 

• Crack extension (Aa) 

• Crack growth rate (da/dN) 

• Number of delay cycles (%) 

The shape of the overload-af fected da/dN versus Aa curve, as well as crack- 
tip residual stresses were also examined. 
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4.1 Influence of Crack Growth Rate Properties on Analytical Predictions 

The FAST-2 model computes AK e ff for each cycle, or block of cycles, 
and then integrates the crack growth rate properties, da/dN as a function 
of AK e ff, to obtain da/dN versus Aa. A procedure which may be used to 
obtain the fatigue crack growth rate properties in terms of AK e ff is de- 
scribed in Reference 7. This procedure involves using the closure model 
to analyze the test specimen from which the crack growth rate data are 
obtained thereby computing AK e ff for each AK and R value. Growth rate 
data at several values of R should in theory collapse into a single curve 
if da/dN is plotted vs AK e ff. 

For crack growth rate data obtained at high-R, little or no crack 
closure should occur; therefore, AK e ff should be nearly equal to AK. Thus, 
the single curve which represents da/dN vs AK e ff for all load ratios should 
then be the same as the high-R da/dN vs AK curve. Thus, if data at high-R 
are available, the procedure described in Reference 7 is not necessary. 

Such data were available from Alcoa from separate production lots 
of 7091-T7E69, as shown in Figure 39. Between Extrusions 1 and 2 (shown 
as dots in Figure 39) there is little variability at R = 0.33. The third 
extrusion from another lot was used at SwRI for the experiments described 
herein. Several crack growth rate measurements at AK ~ 6 MPav^n and R 

= 0.33 were made at SwRI which showed agreement with Extrusions 1 and 

2. However, a fourth extrusion tested by Alcoa at R = 0.33 shows signifi- 
cantly higher growth rates than Extrusions 1 and 2 at R = 0.33 and AK 

< 2 MPat'm. Therefore, there is some uncertainty as to which growth rate 
properties should be used in the model, particularly in the low-growth-rate 
regime, if predictions are to be compared with experiments. 

Figure 39 also shows the range of crack growth rates covered in 
the overload experiments that were analyzed. Very high growth rates are 
obtained during the overload, and very low growth rates are obtained from 
subsequent retardation of the crack growth rates. Another bar on the 
right shows the range of crack growth rates for which constant amplitude 
high-R crack growth rate data were available. Note the unavailability 
of data above 2 x 10"^ m/cycle; this gives rise to uncertainty as to what 
growth rate properties should be used in the model to predict crack exten- 
sion during the overload cycles. 



Crack growth rate (m/cy) 
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Figure 39 shows three curves which were used in the analyses to 
examine the effects of this uncertainty at both very high and very low 
growth rates. Curve C represents a power law fit to the R = 0.8 data 
for Extrusion 2 (shown as x' s in Figure 39). Analyses performed using 
Curve C revealed that the increment of crack extension for the overload 
cycle was much less than that measured in the experiments. The high-AK 
measurements made in the overload experiments indicated that the high 
growth rates were better represented by Curves A or B in Figure 39. 

In the low growth rate regime, Curves A and C are nearly coincident. 
Curve B was intended to estimate the R = 0.8 data for Extrusion 4. For 
Extrusion 2, there exists a consistent offset in AK at each crack growth 
rate between the R = 0.8 and R = 0.33 data. This offset may be thought 
of as the difference between AK and AK e ff. The particular offset at each 
crack growth rate was applied to the R - 0.33 data for Extrusion 4 to 
obtain the portion of Curve B for da/dN < 7 x 10"^ m/cycle. 

The sensitivity of the analysis to differences in constant amplitude 
crack growth rate properties is illustrated in Figures 40 through 43. 
The results in Figure 40 are for Experiment 1 with a constraint factor 
(a) of 1.9. (This constraint factor gave optimum results for this experi- 
ment, as discussed in Section 4.2.) As indicated, the minimum growth 
rate is primarily affected, while the differences in delay cycles is mini- 
mal. Note that the experimental observations indicate the minimum growth 
rate is higher than predicted, and that the measured results recover from 
the minimum rate more gradually than predicted by the model. 

Figure 41 makes the same comparison for Experiment 4 with a = 2.3. 
(This constraint factor gave optimum results for this experiment.) Note 
that because a much greater range in growth rates occurred in Experiment 
4, the differences in the analytical results for different growth rate 
curves are more significant than those of Experiment 1. Curve B gives 
the best prediction of the experimental results. Curve C underestimates 
the maximum growth rates in all cases. Unlike Experiment 1, the initial 
growth rate in Experiment 4 is not well predicted by Curves A or B. 

Figures 42 and 43 show that the initial growth rate before the 
overload in Experiment 3 was also not well predicted by Curves A and B. 
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It is suspected that these initial growth rates may have been influenced 
by transient effects associated with load shedding prior to the overload 
experiments i.e., steady state growth rates may not have been fully achieved 
before the next experiment was conducted. 

Figures 42 and 43 show a comparison of the growth rate effect for 
constraint factors of 1.9 and 2.3, respectively. Note that the sensitivity 
of predictions to growth rate effect is more significant for a - 2.3 (Figure 
43). 

To assess the effect of crack growth rate properties and constraint 
factor, many analyses were performed for each of the three experiments 
using the three growth rate curves described above and constraint factors 
in the range from 1.5 to 3.0. More detailed results from these analyses 
are presented in tabular form in Appendix B. 

The FAST-2 model predicts the general trends in crack growth rate 
response following the load excursions. However, a rigorous assessment 
of the model is limited by its sensitivity to variability in constant 
amplitude growth rate properties. 


4.2 Influence of Constraint Factor on Analytical Predictions 

Analyses were performed for the three crack growth rate curves 
in Figure 39 and for a variety of constraint factors. The constraint 
factor, which ranges from 1.0 for plane stress to 3.0 for plane strain, 
serves to elevate the tensile flow stress, thereby simulating the effect 
of through-thickness constraint on the crack-tip stresses. Options in 
the FAST-2 computer program allow the constraint factor to be fixed at 
a given level (used for most cases) or allowed to vary according to the 
instantaneous ratio of the plastic zone size to the specimen thickness 
(p/t); specifically, for the i t ^ cycle: 


M i 2 

a. = 3.67 - 6.67 — a 7 , 
x t 


( 3 ) 


1.0 < a. ^ 3.0 

i 


where: 



68 


For Experiments 3 and 4, the above relation gave constraint factors 
for the overload cycles of 1-68 and 1.75, respectively. For all other 
baseline cycles, and for the overload cycle of Experiment 1 (OLR* = 2.15), 
the above relation gave a ' constraint factor of 3.0. Fixed constraint 
factors used in the analyses ranged from 1.5 to 2.7. 

Figures 44 and 45 show analytical predictions for Experiment 1 

for various constraint factors and Growth Rate Curves B and C, respectively. 
Best results for Experiment 1 appear to be obtained with a constraint 
factor of between 1.7 and 1.9, although the lower constraint factor predicts 
a minimum growth rate lower than that observed in the experiment. This 
lower a compensates for the fact that the measured minimum growth rate 
persists over a longer increment of crack extension than the predicted 
rates and thereby gives the best prediction of Np. 

Figures 46 and 47 show results at various a*s for Experiment 4 

with Growth Rate Curves B and C, respectively. Because of the large OLR* 
for Experiment 4 and the attendant crack growth rate retardation, for 

a < 2.3, the model generally predicted arrest (i.e., the model was run 

for 90,000 cycles and then aborted before the growth rate recovered). 
Predictions using Curve B and a = 2.1 eventually recovered to give Nj) = 
65 kc. 

The above results for Experiment 4 would seem to be in disagreement 
with the results for Experiment 1, where the optimum fixed a was between 
1.7 and 1.9. The best agreement between measured and predicted results 
for Experiment 4 was obtained with Growth Rate Curve B for the variable-a -- 
for which the constraint factor was 1.68 for the overload and 3.0 for 
the baseline cycling. Because the higher constraint factors gave more 

accurate predictions of the minimum growth rates for Experiment 1, it 
is thought that the use of the low constraint factor for Experiment 1 
acts to artificially depress the minimum growth rate to compensate for 
an inaccurate overprediction of the acceleration, or rate of recovery. 

The results for Experiment 4 for a - 2.5 do not appear to differ 
significantly, especially in Figure 47, indicating the model is more sensi- 
tive to changes in a < 2.5 than to changes in a > 2.5. The variable a 
seems to best represent the rate of recovery throughout the analysis and 
hence the shape of the da/dN-Aa curve for Experiment 4 (Figure 46). 
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Aa From Start of Overload Cycle (mm) 

FIGURE 44. COMPARISON OF EXPERIMENTAL RESULTS AND MODEL PREDICTIONS FOR VARIOUS 
CONSTRAINT FACTORS (a). Uncertainty in experimental results is 
shown as the hatched region. 
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Aa From Start of Overload Cycle (mm) 


FIGURE 45. COMPARISON OF EXPERIMENTAL RESULTS AND MODEL PREDICTIONS FOR 
VARIOUS CONSTRAINT FACTORS (a). Uncertainty in the experi- 
mental results is shown as the hatched region. 










Aa From Beginning of Overload Cycle (mm) 


FIGURE 47. COMPARISON OF EXPERIMENTAL RESULTS AND MODEL PREDICTIONS 
FOR VARIOUS CONSTRAINT FACTORS (ct). Uncertainty in ex- 
perimental results is shown as the hatched region. 
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Figures 48 and 49 show the results of predictions for Experiment 
3. Considering that all conditions for this experiment were the same 
as for Experiment 4, with the exception of the underload , it is seen that 
the underload has a significant effect on diminishing the retardation 
due to the overload, and that the plastic-wake closure model predicts 
this phenomenon. 

Figure 49 shows the results of predictions based on Growth Rate 
Curve C for a from 1.5 to 2.5. Only the lowest a's correctly predict 
the amount of crack extension at recovery, but Np for the low a * s is higher 
than Np observed in the experiment. The optimum a using Curve C seems 
to be about 1.9. The maximum growth rates are significantly underpredicted 
by Curve C. 

Results based on Curve B are shown in Figure 48. The initial growth 
rates are poorly predicted by Curve B. The optimum a for these predictions 
is difficult to judge for a = 1.7 - 2.1; the correct crack extension is 
predicted but the Nj) is overpredicted. For variable a and fixed a = 2.5 ~ 
2.7, the Nj) and da/dN average are well predicted. Note the marked differ- 
ence in the shape of the recovery curve for the variable a . 

4 . 3 Summary 

The optimum constraint factor (a) for each experiment was different 
and dependent on the assumed constant amplitude growth rate properties. 
This is largely because this factor can serve as an adjustment of the 
analysis to compensate for other problems in the prediction, e.g., inaccu- 
rate crack growth rate properties. Overall, best results were obtained 
with variable a ; it seems that using fixed a results in poor prediction 

of the shape of the crack growth recovery curve. 

Initial crack growth rates before the overload were not well pre- 
dicted by the model using a growth rate curve that otherwise gave the 
best post-overload results (i.e., Curve B). Examination of the experimental 
data led to the hypothesis that in some cases the crack growth rate had 
not fully recovered from transients associated with load shedding prior 
to the overload experiments. Future analysis should include simulation 

of the events prior to the overload. 



da/dN (m/cy) 




da/dN (m/cy) 



Aa From Start of Overload Cycle (mm) 


FIGURE 49. COMPARISON OF EXPERIMENTAL RESULTS AND MODEL PREDICTIONS 
FOR VARIOUS CONSTRAINT FACTORS (a). Uncertainty in ex- 
perimental results is shown as the hatched region. 
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The maximum and minimum growth rates predicted in the analysis 
are affected mainly by the da/dN-AK e ff curve, while the amount of crack 
extension during retardation is affected mainly by the constraint factor, 
particularly the constraint ■ factor for the overload cycle which directly 
affects the overload plastic zone size. 
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5.0 DISCUSSION 


It is instructive to further compare the crack opening load values 
obtained locally using the stereoimaging technique with those obtained 
remotely using the elastic compliance technique, as shown in Figure 5. 
The observed differences are at least partially due to the fact that the 
stereoimaging measurements are sensitive to the near-tip response of the 
crack, while the compliance measurements depend only on the global response 
of the cracked specimen and are thus insensitive to local crack-tip 
response. In addition, the latter provides a measure of the average 
through- thickness response of the cracked specimen, while the stereoimaging 
measurements may be dominated by the response of the plane stress region 
at the specimen surface. 

The above implies that the extent of crack closure varies signifi- 
cantly from the surface to the interior of the specimen in proportion 
to the transition in the near-crack-tip stress state from predominantly 
plane stress at the surface to increasingly plane strain toward the inte- 
rior. This view is consistent with the fact that agreement between the 
constant amplitude stereoimaging measurements of AK e ff and those predicted 
using the FAST-2 model can only be obtained when plane stress conditions 
are assumed in the model, Figure 4. The latter is achieved by setting 
the constraint factor to a = 1. Concurrently, the larger AK e f f / AK values 
measured by the compliance technique can be matched by the analytical 
results by using larger a values. 

A three-dimensional view of the crack closure process is consistent 
with recent subsurface measurements in a transparent polymer obtained 
using optical interferometry [8,9] and in a steel obtained using a compli- 
ance technique involving hole drilling along with a special displacement 
transducer [10-12]. These results are also consistent with the measured 
increase in AK e ff which has been observed when the plane stress surface 
regions have been removed by machining subsequent to constant amplitude 
fatigue crack growth [13]. Indirect evidence for this effect is provided 
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by the fact that delay periods following overloads tend to increase as the specimen 
thickness is decreased [14-21]. 

In spite of the above observations, a quantitative understanding 
of the role of through-thickness constraint on crack closure, crack growth, 
and load interaction effects has not yet emerged. Specifically, it is 
not possible to generalize the above results with respect to loading vari- 
ables, particularly applied AK, R, and OLR*. For example, most of the 
results cited above were obtained at relatively high AK values where differ- 
ences in through-thickness constraint are accentuated, thus the applica- 
bility of these results to near-threshold fatigue crack growth is uncertain. 

A study utilizing both remote (average) and local (surface) measurements 
of crack closure in complementary fashion, while systematically varying 
thickness and applied loading variables, is needed to clarify this issue. 

Although substantial variability was observed in the measured AK e ff 
values, as shown in Figure 3, there nevertheless appears to be a systematic 
deviation between the measured and predicted values of AK e ff at both low 
and high values of load ratio, Figure 4. The deviation at low-R values 
may be due to the operation of supplementary mechanisms of crack closure-- 
for example, closure induced by accumulation of oxides [22-25] and asperity 
contact due to surface roughness [26-30]. Of the above candidates, oxide- 
induced crack closure is believed to be most likely in the current system. 

Oxides are known to readily occur within cracks in aluminum alloys which 
have been exposed to laboratory air, either during crack extension [25,31] 
or during subsequent storage [31]. On the other hand, the occurrence 
of crack roughness-induced closure is inconsistent with the relatively 
smooth fracture surface morphology which occurs in the 7091 alloy. 

At high-R, the deviation between measured and predicted values 
would appear to be due to another factor. Specifically, it may be due 
to the fact that the crack-tip element size used in the FAST-2 model is 
too large to accurately represent the localized nature of the crack opening 
process, as shown in Figure 2. Thus, the analytical model, which gives 
the response averaged over the entire element, will predict less crack 
closure than measured very near the actual crack tip. While it is not 
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clear whether or not this factor actually predominates at high-R, this 
can be resolved by a systematic examination of the influence of element 
size on the predicted crack closure over the entire domain of relevant 
loading variables. 

Examination of the crack-tip residual stresses caused by the overload 
and overload/underload cycles indicates that the magnitude of the compres- 
sive residual stresses ahead of the crack-tip is approximately limited 
to the yield stress and is independent of the severity of the overload. 
However, the spacial extent of the compressive residual stress field, 
including both the length and breadth, increases as the severity of the 
overload increases. Furthermore, for the case of the underload, following 
the overload, the extent of the compressive residual stress field is re- 

duced, particularly in width. 

Although additional results are needed to demonstrate that the 
above trend extends to a wider range of loading variables, the current 
results nevertheless demonstrate the importance of the residual stresses 
ahead of the crack tip to the growth of cracks subjected to variable-ampli- 
tude loading. 

Considering the crack-tip strain field in greater detail indicates 
that the distribution of strain within the plastic zone is similar both 

before and after the overload, as shown in Figure 26. Thus, a major effect 
of the overload cycle is to alter the magnitude of the crack-tip strain 

field, Figure 24. Moreover, the crack growth rates were always found 
to correlate with the magnitude of the crack-tip strains -- before, during, 
and after the overload, Figure 38. These findings indicate that the mecha- 
nism of crack extension is the same for both the steady state and overload- 
affected regions of crack growth. 

Thus, the problem of predicting crack growth rates under variable 
amplitude loading appears to be one of determining the correct "driving 

force" for crack growth rather than accounting for differences in the 
material response or intrinsic mechanism of crack extension. The similitude 
which appears to be maintained during variable-amplitude crack growth 
suggests the crack "driving force" can be properly characterized in terms 
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of crack closure concepts through AK e ff, provided that all of the physical 
processes contributing to this parameter are considered. Although the 
latter are generally and qualitatively known to include crack-tip residual 
stress, wake plasticity, oxide formation and accumulation, and asperity 
contact, the challenge which remains is to assess their relative importance 
for various growth rate regimes and formulate quantitative model's capable 
of predicting AK e ff and crack growth rate responses. 
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6.0 CONCLUSIONS AND RECOMMENDATIONS 

1. Under constant amplitude loading, the magnitude of the crack- 
tip opening load and AK e ff--as determined local to the crack 
tip by stereoimaging--depend on the applied values of both 
AK and K^x (or AK and R, or K^x and R) as described by 
the following simple relationship for 7091 aluminum: 

eff _ ^ 5 

AK ’ Kaax 

2. The relation between the extent of the crack which remains 
open and the applied load is highly nonlinear. For example, 
the load required to ,, peel n open the last 20 pm of the crack, 
near its tip, is two to three times higher than that required 
to open the crack to within 200 pm of the tip. Thus, the 
measurement of P 0 or AK e ff is highly dependent on the sensi- 
tivity of the measurement technique. 

3. Local surface measurements of crack- tip opening load using 

the stereoimaging technique give higher values of P Q , and 

correspondingly lower values of AK e ff, than do remote measure- 
ment techniques such as elastic compliance. 

4. The magnitude of the overload/underload cycle controls the 

spacial extent of the compressive residual stress field 
both parallel and perpendicular to the direction of crack 
growth, whereas the magnitude of this compressive stress 

field is limited to the flow stress of the material. These 
results are consistent with the degree of crack growth retar- 
dation following the overload/underload cycle and thereby 

attest to the importance of crack-tip residual stress during 
variable-amplitude fatigue crack growth. 
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5. A correlation exists between the fatigue crack growth rate 

and the crack-tip strain measured before, during, and after 
the overload, thus indicating that the intrinsic mechanism 
of crack growth is strain controlled and is the same for 
constant-amplitude and variable-amplitude fatigue crack 
growth. ' 

6. Similitude in crack- tip plasticity is maintained during 
variable-amplitude fatigue crack growth. This similitude 
is manifested by a uniqueness in (a) strain distribution 
within the crack-tip plastic zone, and (b) the relationship 
between crack-tip strain and crack- tip opening displacement. 

7. It follows from Conclusions 5 and 6 that the problem of 
predicting variable-amplitude fatigue crack growth is one 
of determining the proper "driving force" for crack growth. 
Incorporating crack closure concepts into an effective stress 
intensity factor appears to be a viable approach to character- 
izing the crack "driving force," provided the contributions 
of all relevant physical processes are included. 

8. The FAST-2 analytical model provides reasonable predictions 

of the crack growth rate trends following single overloads 

and overload/underload cycles. However, the sensitivity 
of predictions to the constant-amplitude fatigue crack growth 
rate properties, combined with the lot-to-lot variability 
in these properties, currently limits a rigorous assessment 

of this model. This limitation can be overcome by generating 

constant-amplitude crack growth rate data from the same 

7091 extrusion used for the overload experiments. 

9. The maximum and minimum growth rates predicted by the model 
during the overload-affected crack growth period are sensitive 
to the input growth rate properties, while the extent of 
the overload-affected region is sensitive to the constraint 
factor applied to tensile yielding since the latter affects 
the overload plastic zone size. The shape of the growth 
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rate versus crack extension curve following the overload 
is most accurately predicted by using the variable constraint 
feature of the model. 

10. It is recommended that further assessments of the FAST-2 
analytical model consider the following factors: 

(a) Sensitivity of results to crack-tip element size 

(b) Differences between displacement-based and K-based 
crack opening criterion 

(c) Limitations of applying the Dugdale solutions 
for the center-cracked-tension geometry to other 
geometries*; specifically, by defining the amount 
of crack extension which can be tolerated without 
introducing significant errors due to differences 
in K-gradient, dK/da. 


* The current version of FAST-2 uses the Dugdale solution for a center- 

cracked- tension geometry. In order to properly apply this model to 

the analysis of the single-edge-notched specimens used in this study, 
the K values in the model and the experiment were initially matched 
and the extent of crack extension was limited to that corresponding 
to a single overload event. This restriction precluded analyzing a 
series of overload events to examine the possibility of interactions 
between events. Such analyses would provide useful guidance for the 
planning of experiments since multiple tests, performed on a single 

specimen, are often utilized to maximize the efficiency of data 

generation. 
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APPENDIX A 

DEFINITION OF LOADING VARIABLES AND DELAY CYCLES 
FOR VARIABLE AMPLITUDE FATIGUE CRACK GROWTH 


The quantities related to constant load-amplitude cycling of fatigue 
cracks are given in Fig. A1. Effective AK is computed from values of the load 
P Q required to open the crack down to the tip and the R ratio, according to 
the equation in the figure. 

Several different definitions have evolved to characterize the magnitude 
of the load excursions. Experiments are commonly conducted so that the mean 
stress of the baseline cycling is maintained constant during a given overload 
test by fixing the load ratio (R = = p m i n / p max )- The definition of 
the overload, or underloading ratio can be formulated in terms of either the 
ratio of the maximum values or the ranges of the stress intensity factor for 
the overload to base loading cycles, as shown in Fig. A2, giving 


OLR = %/K,^ = p 0 L /p max 


(A1) 


OLR* = AK ql /AK = 
ULR* = 


(P 0L 


AK 


P • )/(P - P • ) 

r min //vir max r mm y 

P - P 

max UL 

P - P . 

max min 


(A2) 


where P m ^ n and P^x are the minimum and maximum load values in the base cycle, 
and P QL is the maximum load in the overload cycle, and P^ is the minimum load 
in the underload cycle. As indicated in Equations A1 and A2, OLR is refer- 
enced to zero load, while OLR* and ULR* are referenced to the minimum cyclic 
load. Thus, both definitions are identical when P m ^ n = 0 (R = 0). However, 
in general the relationship between OLR and OLR* depends on R as follows: 


OLR - R 


1 - R 


(A3) 


The current study uses OLR* and ULR*, since these definitions are more mean- 
ingful when comparing results at different R values. 
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The various definitions which have been used for the number of delay 
cycles are illustrated in Figure A3. A typical response of crack length ver- 
sus number of cycles and corresponding crack growth rate versus number of 
cycles are shown in Figures 3(a) and 3(b), respectively. The response of the 
crack during a single overload experiment is as follows. Initially, steady 
state crack growth occurs in region a-b at stress intensity factor AK^ imme- 
diately preceding the overload. The overload cycle corresponding to AKq^ may 
then cause a brief accelerated growth period in region b-c, followed by a pre- 
cipitous decrease in growth rate to a minimum growth rate and eventual re- 
covery in region c-e. Steady state growth is re-established in region e-f 
at AK 2 . Generally, AK 2 is nearly equal to AK^ since growth has only occurred 
over a crack length interval which is on the order of the plastic zone size of 
the overload. 

The "delay” caused by the load excursion can be operationally defined in 
terms of a period over which the effective crack growth is zero by construc- 
ting c-d and e-d to give N' D as shown in Figure A3(a). Alternatively, delay 
may be simply defined as the number of cycles over which the crack growth rate 
is less than the preoverload value— this measure corresponds to b-e in Figures 
A3(a) and A3(b) and is labeled Np. 

The above two definitions can differ significantly, particularly for low 
overload ratios where the delay period is relatively small. Np will always be 
greater than N' D since the latter is based on the extrapolation e-d and does 
not include the accelerated growth period b-c. 

The current study uses N D as the measure of delay primarily because it 
is a relatively straightforward operational definition which can be easily ap- 
plied to both experimental data and model predictions by monitoring the post- 
overload crack growth rate. 
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FIGURE Al. DEFINITION OF LOAD AND STRESS INTENSITY FACTOR VALUES 
RELATED TO THE MEASUREMENT OF CRACK OPENING LOAD. 
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FIGURE A2. CHARACTERIZATION OF LOAD HISTORY AND DEFINITIONS 
OF OVERLOAD AND UNDERLOAD MAGNITUDES. 
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APPENDIX B 

ADDITIONAL ANALYTICAL RESULTS OBTAINED USING 
. THE FAST-2 MODEL 


This appendix contains results which supplement those given in Section 
4.0. First, a detailed summary of the comparisons of analytical predictions 
and experimental results is provided in tabular form. Then, the residual 
stresses at the maximum and minimum loads in the cycle are presented at vari- 
ous times before, during and after the overload/underload events. 

B1. Detailed Comparison of Analytical and Experimental Results 

The results for each case studied are presented in Tables B1 through B3. 
Information contained in these tables includes the predicted growth rate and 
effective load ratio (P„/P mov ) prior to the overload (initial) and during the 
overload (maximum) . Also shown for the latter case is the increment of crack 
extension attributable to the overload cycle. The data labeled minimum corre- 
spond to the point of minimum crack growth rate following the overload. The 
minimum growth rate, R e ff minimum, total crack extension from the beginning 
of the overload cycle, and the number of cycles applied between the overload 
and minimum are also shown. Finally, at recovery (i.e., the point when the 
crack growth rate recovers to the initial value), the total crack extension 
and number of cycles (N^) are shown. In order to isolate the error in the 
predicted crack extension due to uncertainty in high growth rates associated 
with the overload cycle, the quantity Aa-Aag^ (i.e., the Aa recovery - Aa max- 
imum) is shown. 

The right-most column in each table reports the experimental data. 
These data are sometimes reported as a range to reflect the uncertainty in the 
experimental measurements. 

Comparing experimental results and analytical predictions reveals that 
the constraint factor and growth rate curve which yields the best estimate of 
Np does not always yield the best estimate of crack extension. Therefore, 
another criterion for comparison is offered in the last row, i.e., the average 
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0.20 

0.010 

10.3 

0.19 

0.011 

11.2 

0.10 

0.011 

12.1 

0.17 

0.012 

14.2 

0.15 

0.011 

15 - 70 
N/A 

0.015 - 0.020 

w da/dN (* 10~ 9 m/cy) 

(1.51)2 

(?.90) 

(5.11) 

(10. 1) 

12.0 

(3.99) 

(5.11) 

4.04 

5. 58 

6.03 

3.63 

(6.39) 

(7.50) 

7.33 

7.90 

0.39 

9.27 

10.0 

1 

N «nff 

(0.7B) 

(0.74) 

(0.60) 

(0.60) 

0.58 

(0.75) 

(0.77) 

0.73 

0.71 

0.69 

0.03 

(0.75) 

(0.72) 

0.73 

0.71 

0.69 

0.67 

0.78 

* A* (»»") 

(0.01?) 

(0.035) 

(0 025) 

(0 019) 

0.020 

(0.036) 

(0.010) 

0.032 

0.030 

0.020 

0.019 

(0.036) 

(0.010) 

0.012 

0.020 

0.020 

0.026 

0.00 

U H (U) 
M 

(0.1) 

(5.2) 

(1.7) 

(0 81) 

0.63 

(19) 

(3.6) 

2.9 

2.2 

1.0 

6.2 

(36) 

(2.7) 

2.1 

1.5 

1.5 

1.0 

4.5 

R 

E Aa (xm) 

('007)3 

(o.m/) 

(0 175) 

(0.120) 

0.117 

(o.ino) 

(0.701) 

0.197 

0.146 

0.146 

0.271 

(0.186) 

(0.700) 

0.109 

0.115 

0.150 

0.143 

0.20 

n « D M 

(>?0.0) 

(15.0) 

(0.4) 

(3 6) 

2.7 

(13.0) 

*9 1) 

0.1 

5.4 

4.3 

19.0 

(10.0) 

(7.1) 

6.5 

4.4 

3.6 

2.5 

0.5 

If Aa-ftam (mm) 

_ 

(0.1(11) 

(0.172) 

(0.126) 

0.114 

(0.170) 

(0.191) 

0.106 

0.135 

0.134 

0.210 

(o.ini) 

(o.ino) 

0, 170 

0.134 

o.no 

0.129 

0.10 0.19 

jj avg. da/dN following 01 
Y (x 10 8 m/cy) 

" 

(i.e» 

(?.o) 

| 

(3.5) 

4.2 

03) 

(2.1) 

2.3 

2.5 

3.1 

l.l 

(i.o) 

(2 5) 

2.7 

3.0 

3.8 

5.1 

1.0 ‘ 2.3 


1. In Ihnse cases wlicre the dashes ate shown. Information on the maximum qrtiwlh rale was not ohtalocd. 

Z. In those cases where results are In parentheses, the number ol elements used In the analysis changed during the confutations and 
Influences the results significantly, confounding comparisons wllh other cases where the number of elements was constant. 

3. If results are ptesented as an Inequality, the computation was stopped at a preset number of cycles (usually 90,000); l.e., the 
model effectively predicts crack arrest. 

4 . Ihe constraint factor was 3.0 (plane strain) throughout the analysis. 



TABLE B2 


EXPERIMENT 3: OLR* = 3.0, ULR* = 2.0, R = 0.5 




GROWTH RATE CURVE 

C 


GROWTH RATE CURVE A ~] 

GROWTH RATE CURV 

E 8 



EXPERIMENT 

I CONSTRAINT FACTOR— 

1.5 

1.7 

1.9 

2.3 

2.5 

2.7 

1.7 

1.9 

2.1 

2.3 

2.5 

2.7 

1.5 

1.7 

1.9 

2.1 

2.3 

2.5 

msm 

Var.2 

I 

N 

I 

da/dN (x 10-7 ra/cy ) 

1.01 

1.22 

1.39 

1.67 

1.75 

1.89 

4.13 

5.01 

5.46 

6.54 

6.93 

6.94 

3.01 

4.13 

5.01 

5.47 

6.54 

6.93 

6.94 

6.94 

1.1-1. 5 

T 

I 

A 

L 

R off 

0.60 

0.58 

0.57 

0.55 

0.54 

0.54 

0.59 

0.57 

0.56 

0.55 

0.55 

0.55 

0.60 

0.59 

0.57 

0.56 

0.55 

0.55 

0.55 

0.55 

0.59 

M 

A 

da/dN (x 10' 5 m/cy) 

.1 

1.47 

1.54 

1.65 

1.67 

1.72 

7.36 

0.07 

0.41 

9.18 

9.47 

9.48 

6.61 

7.36 

0.07 

8.41 

9.18 

9.47 

9.48 

9.48 

3. 0-7.0 

X 

I 

R eff 

- 

0.29 

0.28 

0.27 

0.27 

0.27 

0.29 

0.29 

0.28 

0.28 

0.27 

0.27 

0.30 

0.29 

0.29 

0.28 

0.20 

0.27 

0.27 

0.27 

0.30-0.32 

M 

U 

M 

Aa (mm) 

- 

0.015 

0.015 

0.017 

0.017 

0.017 

0.074 

0.001 

0.084 

0.092 

0.095 

0.095 

0.067 

0.074 

0.081 

0.084 

0.092 

0.095 

0.095 

0.095 

0.03-0.07 

It 

I 

N 

da/dN (x 10‘1° m/cy) 

7.10 

10.9 

17.7 

35.7 

45.5 

56.6 

7.56 

7.04 

7.94 

7.87 

7.97 

10.2 

9.21 

9.65 

9.50 

9.67 

9.92 

9.99 

12.1 

44.1 

15.0-20.0 

R cff 

0.80 

0.78 

0.75 

0.70 

0.60 

0.66 

0.81 

0.81 

0.00 

0.81 

0.81 

0.79 

0.81 

0.81 

0.81 

0.01 

0.81 

0.81 

0.79 

0.72 

0.76 

1 

M 

Aa (nm) 

0.159 

0.105 

0.093 

0.065 

0.064 

0.060 

0.173 

0.159 

0.151 

0.135 

0.131 

0.126 

0.145 

0.173 

0.158 

0.140 

0.139 

0.132 

0.125 

0.231 

0.175 

U 

H 

M (kc) 

12.0 

4.7 

2.7 

0.09 

0.74 

0.57 

9.3 

6.0 

5.0 

3.2 

2.5 

1.6 

5.6 

7.6 

5.5 

3.7 

3.1 

2.2 

1.4 

1.9 

2.6 

R 

Aa (mm) 

— 1 

0.332 

0.242 

0.200 

0. 147 

0.145 

0.150 

0.370 

0.328 

0.282 

0.260 

0.240 

0.220 

0.443 

0.370 

0.325 

0.202 

0.260 

0.239 

0.219 

0.363 

0.29-0.40 

E 

C 

NO (kc) 

19.0 

9.7 

5.2 

2.1 

1.6 

1.3 

16.0 

12.0 

8.8 

7.1 

5.4 

3.7 

19.0 

13.0 

9.9 

7.6 

6.1 

4.7 

3.3 

2.6 

3. 5-5.0 

0 

V 

Aa -Aaql (nm) 

0.317 

0.227 

0.105 

0.130 

0.128 

0.133 

0.296 

0.217 

0.198 

0.168 

0.145 

0.125 

0.376 

0.296 

0.244 

0.198 

0.168 

0.144 

0.124 

0.268 

0.22-0.37 

E 

R 

avg. da/dN following 01 

1.7 

2.3 

3.6 ' 

6.2 

0.0 

10.2 

1.9 

2.1 

2.2 

2.4 

2.7 

3.4 

2.0 

2.3 

2.4 

2.6 

2.8 

3.1 

3.8 

10.3 

4.4-10.6 

Y 

(x 10* 8 m/cy) 







1 





1 












1. In those cases where the flashes are shown. Information on the maximum growth rate was not obtained. 

2. The constraint factor for this case was 3.0 (plain strain) for all cycles except the overload cycle where the 
constraint factor was computed to be 1.75. 











TABLE B3 


EXPERIMENT 4: OLR* = 3.0, R = 0.5 




GROW ill RAIE CURVE C 


GROW HI 


COHSIRAINf T ACTOR 

1.9 

2.3 

2.5 

2.7 

1.9 

2.1 

1 

N 

1 

da /till (x 10*7 m/r.y) 

1.40 

1.67 

1.75 

1.89 

5.01 

5.50 

1 

1 

R cff 

0.57 

0.55 

0.54 

0.54 

0.57 

0.56 

A 








L 








n 

A 

da/dN (x 10' 5 m/cy) 

1.55 

1.65 

1.69 

1.74 

0.16 

8.50 

X 

1 

R eff 

0.28 

0.27 

0.27 

0.27 

0.29 

0.20 

M 

U 

Aa (mm) 

0.016 

0.017 

0.017 

0.017 

0.082 

0.085 

M 








M 

1 

da/dN (x 10-10 m/cy) 

<3.8/1 

5.07 

0.75 

12.9 

<0.52 

<1.35 

N 

R eff 

>.92 

0.91 

0.90 

0.09 

>0.96 

>0.95 

1 

M 

Aa (nm) 

>0.042 

0.015 

0.013 

0.041 

>0.115 

>0.125 

U 

M 

N (kc) 

>20.0 

21.0 

14.0 

9.4 

>90.0 

>90.0 

R 

E 

Aa (nm) 

>0.012 

0.100 

0.143 

0.147 

>0.115 

>0.125 

C 

0 

H 0 

>20.0 

12.0 

25.0 

24.0 

>90.0 

>90.0 

V 

Aa-AaoL (nw) 

- 

0.171 

0.126 

0.130 

- 

- 

E 

D 

avg. da/dN following OL 

- 

0.41 

0.50 

0.51 

- 

- 

K 

Y 

(x I0* 8 m/cy) 








ATE CURVE A 1 

GR0W1II RAf 

E CURVE D 



2.3 

2.5 

2.7 

1.9 

2.1 

2.3 

2.5 

2.7 

Yar.* 

6.53 

6.90 

7.00 

5.05 

5.51 

6.53 

6.97 

6.90 

6.98 

0.55 

0.55 

0.55 

0.57 

0.56 

0.55 

0.55 

0.55 

0.55 

9.21 

9.57 

9.58 

0.16 

8.50 

9.21 

9.57 

9.58 

9.57 

0.28 

0.2/ 

0.27 

0.29 

0.28 

0.20 

0.27 

0.27 

0.27 

0.092 

0.096 

0.096 

0.002 

0.089 

0.092 

0.096 

0.096 

0.096 

7.68 

24.9 

37.7 

<0.71 

4.42 

25.7 

44.9 

61.4 

17.66 

0.92 

0.90 

0.88 

>0.97 

0.96 

0.93 

0.90 

0.88 

0.94 

0.126 

0.126 

0.123 

>0.121 

0.123 

0.125 

0.127 

0.123 

0.169 

17.0 

6.9 

4.5 

>90.0 

24.0 

6.9 

4.5 

3.1 

14.0 

0.336 

0.310 

0.277 

'0.121 

0.319 

0.336 

0.310 

0.276 

0.449 

49.0 

20.0 

13.0 

>90.0 

65.0 

22.0 

13.0 

9.3 

45.0 

0.241 

0.214 

0.181 

- 

0.230 

0.214 

0.214 

0.100 

0.353 

0.50 

l.l 

1.4 


0.35 

1.1 

1.6 

1.9 

0.78 


1. |f results are presented as an Inequality, the computation was stopped at a preset number of cycles (usually 90.000); l.e.. the 
model effectively predicts crack arrest. 

2. the constraint factor for this case was 3.0 (plain strain) for all cycles except the overload cycle where the constraint factor 
was computed to be 1.60. 


EXPERIMENT 


1.5 - 1.7 
0.70 

5.0 - 0.0 
0.25 

0.05 - o.on 
i.o - lon.o 

0.96 

0.10 - 0.14 

3.0 - 7.0 


>0.3 

>16.0 

0.22 - 0.35 
1.1 - 2.1 
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crack growth rate following the overload which is the quantity Aa - AaQ^ di- 
vided by Np. 

B2. Residual Stresses Predicted Before, During and After the Overload/ 

Underload Event 

♦ 

B2.1 Experimental 1: OLR* = 2.15, R = 0. 1 

Stresses computed by the model are shown before, during and after the 
overload cycle in Fig. B1. These results are independent of the crack growth 
rate curve used (except for the number of cycles), but depend on aK and the 
overload ratio used, as well as the constraint factor (a = 1.9). 

The magnitude of stress at maximum load (solid line) is equal to the 
constraint factor times the flow stress, while the magnitude of residual 
stress at minimum load (broken line) is limited to the unconstrained flow 
stress. Elastic stresses are not plotted, but return to average applied 
stress in a manner proportional to the square root of the distance from the 
crack tip. The stair step appearance of the stress is due to the division of 
distance into discrete elements. Stress is shown as constant for the width of 
each element. 

Prior to the overload, Fig. Bla, the extent of the residual stresses at 
minimum load is small compared to that on the overload cycle, Fig. Bib. Ten 
cycles after overload, Fig. Blc, the residual stress fields at both maximum 
and minimum load have been modified considerably, as compared to the zone 
prior to overload, in that there is now a region of minimum tensile stress 
within the plastic zone at maximum load. Also, the extent of residual stress 
in the crack wake at minimum load has been diminished, compared to that prior 
to the overload. 

The crack is seen to have grown into the plastic zone of the overload 
3513 cycles later, Fig. Bid. This is near the point of minimum growth rate. 
Stress at both maximum and minimum load has been considerably altered, and 
there is again residual stress in the crack wake at minimum load. Residual 
stress ahead of the crack at minimum load is seen to be tensile, as compared 
to compressive before the overload. 



Stress/Flow Stress 




Olstance from Crack Tip (mm) Distance from Crack Tip (mm) 


d) N = 3513, a = 4.18 mm e) N = 5035, a - 4.20 mm 

(minimum growth rate) 



Distance from Crack Tip (mm) 

f) N = 7329, a = 4.32 mm 
(recovery) 


Figure Bl. Predicted crack-tip stress distributions and crack closure levels for cycles* (N) before, during 
(N=0), and after the overload in Experiment 1 at AK = 7 MPa/in, R “ 0.1, OLR* “ 2.15* Curves 

correspond to maximum ( ) and minimum ( ) applied load. Predictions for a *= 1.9 for all 

cycles . 
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At 5035 cycles, the condition is altered, Fig. Ble, but not appreciably, 
over that at 3513 cycles, but the crack has grown only about 0.1 mm. However, 
after 7329 cycles, Fig. Blf, considerable crack growth has occurred, and the 
residual stress patterns at both maximum and minimum loads have come to look 
more like those before application of the overload. In addition, there is now 
a small region of compressive residual stress in the crack wake just ahead of 
the spot where the overload was applied. 

B2.2 Experiment 3: OLR* = 3.0. OLR* = 2.0. R = 0.5 

Comparison of Fig. B2a, one cycle prior to overload, with the same 
condition at R = 0.1, Fig. Bla, shows that there is significantly less 
residual stress in the crack wake at minimum load for this high R. Note that 
a constraint factor of 3 is being used in this case, as compared to the 1.9 
used for Experiment 1 . 

The overload followed by underload results in a large zone of residual 
stress at minimum load which extends both ahead and behind the crack tip. The 
residual stress zone is much larger for this experiment than for Experiment 1 , 
even though the maximum load is about the same, because the crack is over 
twice as long. Note that the constraint factor has also changed for this 
cycle from 3 to 1 (see Fig. B2b). 

Two cycles after the overload/underload, the crack has advanced approxi- 
mately 0.1 mm, Fig. B2c. Although the residual stresses are similar to those 
of Experiment 1, Fig. Bic, their extent is considerably greater. 

After 19^0 cycles, Fig. B2d, near the point of minimum growth rate, a 
considerable wake of residual stress has developed and only a small region of 
compressive residual stresses has developed at the crack tip at minimum load, 
in clear contrast with the result of Fig. Bid. 

After 2600 cycles, the growth rate has nearly recovered to that prior to 
the overload and the residual stress field is also similar, Fig. B2f. 

B2.3 Experiment 4; 0LR« = 3.0. R = 0.5 

Computed residual stress distributions for this case are shown in Fig. 
B3. Prior to application of the overload, the residual stress distribution is 



Stresss/Flow Stress , Stress/Flow Stress 


a) N = -1, a = 10.44 mm 


b) N = 0, a = 10.44 mm 


c) N = 7 ,'a = 10.54 mm 



Distance from Crack Tip (nm) 

d) N = 1940, a = 10.68 mm 
(minimum growth rate) 


Distance from Crack Tip (mm) 

e) N = 2263, a = 10.70 mm 



Distance from Crack Tip (mm) 

f) N = 2601, a = 10.71 mm 
(recovery) 


Figure B2 . Predicted crack-tip stress distributions and crack closure levels for cycles (N) before, during 

(N=0), and after the overload and underload in Experiment 3 at AK 7 MPa/m, R 0,5 » ’ 

ULR* = 2.0. Curves correspond to maximum ( ) and minimum ( ) applied load. Pre ictions or 

a = 3, except during the overload where a = 1.75. 


Stress/Flow Stress Stress/Flow Stress 


-.08 -.04 0 


a) N = -1, a = 8.95 


b) N = 0, a = 8.95 


.4 .5 


Distance from Crack Tip (mm) 

d) N * 14,607, a = 9.12 
(minimum growth rate) 


"- 1 o .1 .2 .3 

Distance from Crack Tip (mm) 

e) N = 31,378, a = 9.16 


-.2 -.1 


Distance from Crack Tip (mm) 

f) N = 45,286, a = 9.31 
(recovery) 


Predicted crack-tip stress distributions and crack closure levels for cycles (N) before, during 
(N=0), and after the overload in Experiment 4 at AK = 7 MPav^n, R = 0.5, OLR* = 3.0. Curves cor- 
respond to maximum ( ) and minimum ( ) applied load. Predictions are for a = 3, except 

during the overload where a = 1.68. 
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similar to those for Experiment 3, which was at the same R. Residual stress 
after the overload cycle is also similar to that of Experiment 3, Fig. B2b, 
but the extent of the stresses is greater and there is no residual stress in 
the wake as there was when the overload was followed by an underload. 

Three cycles after the overload, the residual stress distribution is as 
shown in Fig. B3c, which is considerably different than that shown in Fig. 
B2c, in that the residual stress is always negative ahead of the crack tip, 
rather than having a positive component. 

After 14,600 cycles, the crack has grown well into the zone of residual 
stress and residual stress in the wake has developed, Fig. B3d. This is near 
the point of minimum growth rate. At 31,378 cycles after the overload, Fig. 
B3e, the max imum stress has increased considerably, compared to that in Fig. 
B3d. At 45,286 cycles, the growth rate has nearly returned to that previous 
to the overload, and so has the stress distribution, Fig. B3f. 
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